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Abstract 

We study the Higgs boson pair production through e + e~ collision in the noncommutative (NC) 
extension of the standard model using the Seiberg-Witten maps of this to the first order of the 
noncommutative parameter Q^ u . This process is forbidden in the standard model with background 
space-time being commutative. We find that the cross section of the pair production of Higgs boson 
(of intermediate and heavy mass) at the future Linear Collider(LC) can be quite significant for 
the NC scale A lying in the range 0.5 — 1.0 TeV. Finally, using the direct experimental(LEP II, 
Tevatron and global electro-weak fit) bound on Higgs mass, we obtain bounds on the NC scale as 
665 GeV < A < 998 GeV. 



*Author(corresponding) : pdas@bits-goa.ac.in, pdasMaparna@gmail.com 
^ abhishodh@gmail.com 
1 anupam . mitr a@ gmail . com 



1 



I. INTRODUCTION 



Inspite of its enormous experimental success the standard model(SM) of particle physics 
still awaits the discovery of the Higgs boson. After the Large Electron Positron (LEP) 
collider has set a lower limit of about 114.4 GeV on its mass [l|, the responsibility of finding 
the Higgs now rests mostly on the Large Hadron Collider (LHC) at CERN which has started 
its operation. At the same time, puzzles such as the naturalness problem make a strong case 
for physics beyond the standard model(SM), just around or above the mass scale where the 
Higgs boson is expected to be found. It is therefore of supreme interest to see if the collider 
signals of the Higgs boson contain some imprint of new physics. This necessitates detailed 
quantitaive exploration of a variety of phenomenona linked to the production and decays of 
the Higgs. 

In this paper, we have studied pair production of the Higgs boson in the intermediate 
and heavy mass range at the Linear Collider(LC) as a possible channel for uncovering new 
physics effects. In particular, we show that such pair production which is forbidden in the 
SM with commutative space-time (we will call this as CSM in abbreviation) receives a large 
contribution in the noncommutative(NC) extension of the SM. 

As has been mentioned above, the large hierarchy between the electroweak scale Mw 
and the Planck scale Mpi is somewhat puzzling. Though theories like supersymmetry and 
technicolour, each with its own phenomenological implications and constraints, have been 
proposed as a resolution of this problem, the idea of extra spatial dimensions with the scale 
of gravity being as low as TeV, has drawn a lot of interest among the physics community 
recently |2|. In some brane- world models jsj where this TeV scale gravity is realised, one can 
principly expect to see some stringy effects in the upcoming TeV colliders and in addition 
the signature of space-time noncommutavity. Interests in the noncommutative(NC) field 
theory arose from the pioneering work by Snyder (4j and has been revived recently due to 
developments connected to string theories in which the noncommutativity of space-time is an 
important characteristic of D-Brane dynamics at low energy limit js-7|. Although Douglas et 
al. 6| in their pioneering work have shown that noncommutative field theory is a well-defined 
quantum field theory, the question that remains is whether the string theory prediction and 
the noncommutative effect can be seen at the energy scale attainable in present or near future 
experiments instead of the 4-d Planck scale M v [. A notable work by Witten et al. (9] suggests 
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that one can see some stringy effects by lowering down the threshold value of commutativity 
to TeV, a scale which is not so far from present or future collider scale. 

What is space-time noncommutavity? It means space and time no longer commute with 
each other and as a result one cannot measure the space and time coordinates simultaneously 
with the same accuracy. Writing the space-time coordinates as operators we find 

[x lt ,x v \ = ie lu , (i) 

where the matrix Q^ u is real and antisymmetric. The NC parameter Q^ u has dimension of 
area and reflects the extent to which the space-time coordinates are noncommutative i.e. 
fuzzy. Furthermore, introducing a NC scale A, we rewrite Eq. [1] as 

[X„X V ] = ^ (2) 

where 0^(= c^/A) and has the same properties as Q^ u . To study an ordinary field 
theory in such a noncommutative fuzzy space, one replaces all ordinary products among the 
field variables with Moyal-Weyl(MW) [10] * products defined by 

(f*g)(%) = exp (-Q^d^dy^ f{x)g{y)\ y=x . (3) 

Using this we can get the NCQED Lagrangian as 

C = * - (Drf) * 7 "V) - mj> * V - \f,u * F^, (4) 

which are invariant under the following transformations 

tp(x,Q) ^tfj'(x,&) = U*il){x,Q), (5) 
A„(x, 6) -> A' Ax, &) = U ★ A^x, 6) * U~ l + -U* d^U~\ (6) 



where U = (e iA )*. In the NCQED lagrangian (EqgJ) D^ip = d^ip - ieA^ * {D^ip) = 
+ ieip * Ap, = dpAv - d v Afj, - ie^A^ ★ A. - A v ★ A^). 
The alternative is the Seiberg-Witten(SW)p-l8( approach in which both the gauge pa- 
rameter A and the gauge field A^ is expanded as 

A a (x, 9) = a(x) + 6^A«(x; a) + Q^Q^A^x; «) + ••• (7) 
A p (x, 6) = A p {x) + &»A$ p (x) + Q^Q^A^x) + ■■■ (8) 
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and when the field theory is expanded in terms of this power series Eq. (J7J) one ends up with 
an infinite tower of higher dimensional operators which renders the theory nonrenormaliz- 
able. However, the advantage is that this construction can be applied to any gauge theory 
with arbitrary matter representation. In the WM approach the group closure property is 
only found to hold for the U(N) gauge theories and the matter content are found to be in 



the (anti)-fundamental and adjoint representations. Using the SW-map, Calmet et al. 11] 
first constructed a model with noncommutative gauge invariance which was close to the 
usual commuting Standard Model(CSM) and is known as the minimal NCSM(mNCSM) in 
which they listed several Feynman rules comprising NC interaction. Intense phenomenolog- 
ical searches [3| have been made to unravel several interesting features of this mNCSM. 
Hewett et a/.explored several processes e.g. e + e~ — > e + e~ (Bhabha), e~e~ — > e~e~ (Moller), 
e~7 — > e~7, e + e~ — > 77 (pair annihilation), 77 — > e + e~ and 77 — > 77 in context of NCQED. 
Recently, in a work [13( one of us has investigated the impact of Z and photon exchange in 
the Bhabha and the Moller scattering and three of us have re por ted the impact of space-time 
noncommutativity in the the muon pair production at LC 14J. Now in a generic NCQED 
the triple photon vertex arises to order 0(Q), which however is absent in this minimal 
mNCSM. Another formulation of the NCSM came in forefront through the pioneering work 



by Melic et a/.[15j where such a triple neutral gauge boson coupling jl?j] appears naturally 
in the gauge sector. We will call this the non-minimal version of NCSM or simply NCSM. 
In the present work we will confine ourselves within this non-minimal version of the NCSM 
and use the Feynman rules for interactions given in Melic et a/.[lj]]. 

In Sec. II we present the cross section of e + e~ — > Z — > HH in the NCSM, a process 
which is forbidden in the CSM. The detailed numerical analysis of the pair production cross 
section, angular distribution and the prospects of TeV scale noncommutative geometry are 
presented in Sec. III. Finally, we summarize our results in Sec. IV. 

II. HIGGS PAIR PRODUCTION AT FUTURE LC 

As mentioned before the pair production of a neutral Higgs boson from e + — e~ annihi- 
lation in the CSM is forbidden(at the tree level). So any excess in the predicted event rate 
may be interpreted as the signature of new physics. Supersymmetry and the brane world 
gravity are front runners ( see 16( and references therein). Here we explore the potential 
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of the NCSM. In the NCSM the pair production of a neutral Higgs boson proceeds through 
e + e~ — >■ Z —7- H H via the s channel exchange of Z boson. The corresponding Feynman 
diagram is shown in Fig. [TJ The scattering amplitude using Feynman rules to order for 




FIG. 1: Feynman diagrams for e + e 



HH in the NCSM. 



each vertex can be written as 
%A 



vra M 2 H 



sin 2 (2^)(s — m| + iT z m z ) 



X(j9 3 6)' 



1 + 2^2©Pi) 



(9) 



where s = (pi +P2) 2 , o. = e 2 /47r and 9w is the Weinberg angle. Mh is the Higgs mass and 
rriz and V z are the mass and decay width of the Z boson. The Feynman rules required for 
the above result are listed in Appendix A. The spin averaged squared-amplitude is: 

_ \A\ 2 _ n 2 a 2 M% [l + (4sin 2 fl w -l) 2 l 



\A\ 



sm\26 w ) [(a - m|) 2 + T 



l m V\ 



[2(p 3 &Pi)(p 3 &P2) - (pi.p 2 )[(p 3 e).(p 3 e)]] 



(10) 



Explicit expressions for various dot product terms of the above equation are listed in Ap- 
pendix B. The differential cross-section can be written as 



da 



1 A^M^M 2 )^ 



■\A\< 



dQ 64tt 2 s s 
where a = a(^, A, 8, <p) and A is the Kallen function defined as X(x, y, z) 
2xy — 2yz — 2zx. From [TT1 we can obtain a, da jd cos 6 and da/d(f) as: 



(11) 



x 2 + y 2 + z 2 



a 



da 



dcos9 
da 



1 r2n da 

d(cos6) / d(f)— 
-1 Jo dil 

^ ,da 
d\l 



(12) 
(13) 
(14) 
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III. NUMERICAL ANALYSIS 



In this section, we analyze the total cross section and angular distributions of the cross 
section of the neutral Higgs pair production. Before making a detailed analysis, let us 
make some general remarks regarding the observation of non commutative effects. Since we 
assume = (cQj,Cy) = 6ijkX k ), where & = (E)i and Xk = {B)k are constant vectors in 
a frame that is stationary with respect to fixed stars, the vectors (E)i and (B)k point in fixed 
direction which are the same in all frames of reference. However, as the earth rotates around 
its axis and revolves around the Sun, the direction of E and B will change continuously with 
time dependence which is a function of the coordinates of the laboratory. The observables 
that are measured will thus show a characteristic time dependence. It is important to be able 
to measure this time dependence to verify such non commutative theories. In our analysis, 
we have assumed the vectors E — -^=(i+j + k) and B = -^={i + j + k) i.e. they behave like 
constant vectors. This can be true only at some instant time at most. 

A. Pair production cross section in the NCSM 

As we have mentioned earlier the Higgs pair production is forbidden in the CSM. Any 
signature of such an event will correspond to new physics and NCSM is one of the promising 
candidates among such a class of new physics models. In Fig. (TSJ) we have plotted the total 
cross section er(e~e + — >• Z — >■ H H) as a function of the Higgs mass m^(GeV). The machine 
energy is fixed at E com (= y/s) = 500 GeV and 1000 GeV as shown above. In each figure while 
going from the top to the bottom, the curves correspond to A = 500, 600, 700, 800, 900 
and 1000 GeV, respectively. Note that the pair production cross sections are maximum at 
Tdn = 220(437) corresponding to the machine energy y/s = 500(1000) GeV. Assuming an 
integrated luminsoity of the futute LC about £ = 500 we have predicted the number 

of events in the NCSM. The results, the number of events N (yr _1 ) as a function of A 
corresponding to the machine energy ^/s = 500 GeV and 1000 GeV are presented in Table 
1. Note that if we increase A from 500 GeV to 1000 GeV, the number of events(NC signals) 
A(yr _1 ) decreases from 34(128) per year to 2(8) per year. So a maximum of 34 and 128 
events (NC signal) per year are expected to be observed at the Future LC. 
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FIG. 2: The cross section a(e~e + — > Z 
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HH) (fb) as a function of the Higgs mass m# is shown. 
On the l.h.s(r.h.s) figure center-of-mass energy is fixed at ^fs = 500(1000) GeV. While going from 
the top to bottom the NC scale A increases from 500 GeV to 1000 GeV in steps of 100 GeV. 
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128 
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0.0664 
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500 
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500 


3 
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0.0243 


500 


12 
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0.0042 


500 


2 




1000 


0.0159 


500 
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Table 1: The yearly number of events(NC signals) with the increase in the NC scale A is 

shown. The integrated luminosity of the LC is assumed to be C = 500 fb^ 1 . . 

These are to be compared with the zero event prediction in the CSM where this process is 

forbidden. 
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B. Angular distribution of muon pair production e e + — > Z — > HH in the NCSM 

Azimuthal distribution: The angular distribution of the final state particles is a useful tool 
to understand the nature of new physics. Next we will see how the azimuthal distribution can 
be used to seperate out the noncommutative geometry(and thus the NCSM) from the other 
type of new physics models e.g supersymmetry, brane world gravity, unparticle scenario, 
little Higgs models etc. In Fig. [3] we show 4r as a function of the azimuthal angle <p. For the 




<|> (radian) 4> (radian) 

FIG. 3: The ^{fb/rad) distribution function of <^(in rad) is shown. The lowest horizontal 

curve at zero is due to the CSM, whereas the plots above the horizontal one, as we move up 

correspond to A = 1.0,0.9,0.8,0.7,0.6 and 0.5 TeV, respectively in the NCSM. 

angular analysis study, we consider two cases corresponding to the c.o.m energies -y/i = 500 
GeV and 1000 GeV, respectively. Corresponding to a particular c.o.m energy we vary the 
NC scale A from 0.5 TeV to 1.0 TeV. The distribution da /dip is completely flat at zero in 
the CSM, since the process is forbidden there. Other plots above the horizontal zero line, as 
we move up correspond to A = 1.0, 0.9, 0.8, 0.7, 0.6 and 0.5 TeV, respectively in the NCSM. 
The departure from the flat behaviour is due to P3<9pi, P3OP2 and (p30).(p30) terms in 
Eqns. [10] that bring in the <fi dependence which is observed in Fig. |3j Interestingly, in each 
of the two figures the curves shows maxima and minima. The maxima arises at <fi = 37r/4 
rad, whereas the minima is found to be located at <fi = 77r/4 rad. Note that if we set A = 00, 
the lowest (zero) horizontal curve in the CSM is recovered. 

Note that such an azimuthal distribution clearly reflects the exclusive nature of the non- 
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commutative geometry which is rarely to be found in other classes of new physics models. 

Polar distribution: To extract the exclusive nature of noncommutative geometry, like 
azimuthal distribution, the polar distribution might also be useful. In Fig. HJ -^-^ is plotted 



as a function of cosO. Note that the distribution around the cosO = line of Fig. H] is 

E com = 0.5TeV 

0.02 i , 1 1 1 0.006 



E com =1.0TeV 
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o 
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e 

-a 



-0.5 




cosO 
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FIG. 4: The 3555 (/&) distribution ELS cl function of cosO (in rad) is shown. The lowest horizontal 
curves(at zero value) in both figures are due to the CSM, whereas the plots above the horizontal 
one, as we move up correspond to A = 1.0, 0.9, 0.8, 0.7, 0.6 and 0.5 TeV, respectively in the NCSM. 

completely symmetric. The lowermost plot (zero horizontal curve) in Fig. H] corresponds to 
the polar distribution in the CSM (reflecting the fact that the process is forbidden in the 
CSM) and the plots, as we move up correspond to A = 1.0,0.9,0.8,0.7,0.6 and 0.5 TeV, 
respectively, gets more and more curved, which is a unique features of the NCSM. The 
uppermost curve in the figure corresponding to A = 0.5 TeV exhibits maximal deviation 
from the lowermost CSM curve (obtained by setting A — > 00). 



C. Constraining the NC scale A using experimental bound on Higgs mass ran 

We found in an earlier subsection that the pair production cross section is maximum at 
ran — 220 GeV and 437 GeV corresponding to the c.o.m energy yfs = 500 GeV and 1000 
GeV. The corresponding NC scale is A = 500. In Fig. 5 we make some contour plots in the 
vtlu — A plane corresponding to event production rate per year A^(yr _1 ) with 

• Scenario I : 2 < N < 34 and the c.o.m energy y/s = 500 GeV, 



• Scenario II: 8 < AT < 128 and the c.o.m energy y/s = 1000 GeV. 



E com =0.5 TeV, 2 < events per year <34 
E„ m =1 .0 TeV, 8 < events per year <128 tmmm 
m H = 114.4 GeV - 

m H = 200 GeV 

m H = 158 GeV 

m H = 175GeV 



500 




A(GeV) 

FIG. 5: The contour plot in the ran — A plane obtained by setting 2 < N{yr^ 1 ) < 34(lower) and 
8 < Nd/r^ 1 ) < 128(upper) corresponding to \fs = 500 GeV and 1000 GeV, respectively. The 
lowermost and uppermost horizontal lines correspond to the lower and upper bound on the Higgs 
mass mn which follows from the LEP II direct search of Higgs boson and the global electro-weak 
fit. The combined CDF and DO data at Tevatron excludes 158 GeV < run < 175 GeV. Using 
these we obtain bounds on the bound NC scale A. (See the text for further details). 

The following results are in the order: 

• The direct search of Higgs boson at LEP II gives a lower bound on Higgs mass m# 
as 114.4 GeV. Incorporating this in Fig (jHJ), one finds : (i)A > 666 GeV in scenario I 
and (ii) A > 347 GeV in scenario II. 

• The combined CDF and DO data at Tevatron, Fermilab excludes m# lying between 
158GeV < m H < 175 GeV. Translating this in Fig ©, one finds A < 872 GeV and 
A > 938 GeV in scenario I and A < 468 GeV and A > 515 GeV in scenario II. 

• The Global Electro- weak fit suggests ra# < 200 GeV. Translating this in Fig (jSJ), one 
finds A < 999 GeV in scenario I and A < 581 GeV in scenario II. 

Altogether, the direct searches of Higgs at LEP II and Tevatron and global electroweak 
fit to uih gives rise the following bound on m#: 
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• 665 GeV < A < 872 GeV and 938 GeV < A < 998 GeV in scenario I. 



• 347 GeV < A < 468 GeV and 515 GeV < A < 581 GeV in scenario II. 



IV. CONCLUSION 



The idea of the TeV scale space-time noncommutativity draws a lot of attention in the 
physics community. We explored the impact of spacetime noncommutivity in the Higgs 
pair production e + e~ — > Z — > HH within the NCSM. The plots showing the cross section 
(j(e + e~ — > Z — > HH) as a function of Higgs mass mu at a fixed machine energy E com , sug- 
gests that at a particular m H the pair production cross section is maximum. Corresponding 
to the macine energy y/s = 500 GeV and 1000 GeV, the peak(where the cross section and 
thus the event rate is maximum) is located at m H = 220 GeV and m H = 437 GeV. The 
maximum number of events corresponding to machine energy y/s = 500 GeV and 1000 GeV 
turn out to be 34 (yr^ 1 ) and 128 (yr~ l ), respectively with no events of Higgs pair produc- 
tion in the CSM. Another interesting feature of our study is the angular distribution. The 
azimuthal distribution da/dcj) which is supposed to be completely symmetric (zero in the 
CSM), deviates substantially from its flat (zero) behaviour in the NCSM. The deviation sees 
an enhancement with the lowering of NC scale A. Note that such a non-trivial azimuthal 
distribution is a unique feature of the noncommutative standard model and it is rarely found 
in any class of new physics models(e.g.supersymmetry, Randall-Sundrum model etc) which 
supports neutral Higgs pair production. Finally, we study da/dcosd as a function of cosd: 
zero in the CSM deviates from it's CSM behaviour and develops an asymmetry around the 
cosd = curve in the NCSM. There is a clear departure of the distribution obtained in 
the NCSM from the one obtained in the CSM. We make the contour plots in the plane of 
m H — A corresponding to the event rate 2 < N(yr~ 1 ) < 34 and 8 < N(yr~ v ) < 128. Finally 
the direct searches of Higgs at LEP II and Tevatron and global electroweak fit to m H gives 
rise the following bound on m H : (i) 665 GeV < A < 872 GeV and 938 GeV < A < 998 GeV 
in scenario I and (ii) 347 GeV < A < 468 GeV and 515 GeV < A < 581 GeV in scenario 
II. So, the noncommutative geometry is found to be quite rich in terms of its phenomeno- 
logical implications and it is worthwhile to explore several other processes which might be 
interesting and potentially relevant for the Linear Collider experiments. 
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Appendix A: Feynman rules to order 0(0) 



The Feynman rule for the vertex f(pi n ) — f(p ou t) — Z{k) is [15 ] 



sin29 w 



2sin29 w 



(Al) 



The momentum conservation reads as pi n + k = p ou t- Similarly, the Feynman rule for the 
interaction vertex vertex H(p) — H(p) — Z(k) is: 

g mi (pev 



4cos# 



ir 



(A2) 



Here (pQ)^ = p^Q", T A = (c e v ± c e A j 5 ) . Also p out Qpin = pLt^^Pin = ~Pin®P, 



out- 



Appendix B: Momentum prescriptions and dot products 

Working in the center of momentum frame and ignoring the electron mass, we can specify 
the 4 momenta of the particles as follows: 

» = (f ) (Bl) 



p 3 = (~2~> k' sin9 coscj), k' sinOsmcf), k'cosOj (B3) 
P4 = — k' sin 9 cos (p, — £/ sin # sin 0, — k' cos9 j , (B4) 



2 V s 

where 9 is the scattering angle made by the 3-momentum vector p% of H{p^) with the +ve Z 
axis and is the azimuthal angle. We note that the antisymmetric 0^ has 6 independent 
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components corresponding to c^ u = (co^Qj) with i,j = 1,2,3. Assuming all of them are 
non-vanishing they can be written in the form 



CQi 



A 2 ' 
A 2 



(B5) 
(B6) 



The antisymmetric 0^ is analogous to the field tensor where and Xi are like the 
components of the Electric and Magnetic Field vectors. Setting £j = (E)i = % — 1,2,3 
and Xi = (B)i = 4g, i — 1, 2, 3( noting the fact that Xi — ~X\ 6 — ~C an d 6^ — \$l an d 



XiX j = H> we find 



P1-P2 

p 3 Qpi 
P3OP2 

(p 3 ©)o 

(p 3 e)i 
(p 3 ©) 2 
(p 3 ©) 3 



2' 



2v/3A 2 ' 

— s 
4v/3A 2 

— s 
4^3A 2 



/ 4M 2 

1 -^(2sin#sin0 + cos0) - 1 



/ 4M 2 

1 -^(2sin#cos0 + cos#) + 1 



2v/3A 2 
2V3A 2 



AM 



— (sin 9 cos + sin 9 sin + cos 9) 



AM 



— (cos 9 — sin sin 0) 



4M 



— (sin9 cos — cos 0) 



2^3A 2 



1 + \ 1 



am; 



— (sin9 cos — sin sin . 
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(B8) 

(B9) 
(BIO) 
(Bll) 
(B12) 
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